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In this work, we revisit the electron transfer rate theory, with particular interests in the distinct 
quantum solvation effect, and the characterizations of adiabatic/nonadiabatic and Markovian/non- 
Markovian rate processes. We first present a full account for the quantum solvation effect on 
the electron transfer in Debye solvents, addressed previously in J. Theore. & Comput. Chem. 5, 
685 (2006). Distinct reaction mechanisms, including the quantum solvation-induced transitions 
from barrier-crossing to tunneling, and from barrierless to quantum barrier-crossing rate processes, 
are shown in the fast modulation or low viscosity regime. This regime is also found in favor of 
nonadiabatic rate processes. We further propose to use Kubo's motional narrowing line shape 
function to describe the Markovian character of the reaction. It is found that a non-Markovian 
rate process is most likely to occur in a symmetric system in the fast modulation regime, where the 
electron transfer is dominant by tunneling due to the Fermi resonance. 



I. INTRODUCTION 

A. Prelude 

Electron transfer (ET) is the simplest 
tion system but plays a pivotal role in 
chemical and biological processes. The 
of ET research has grown enormously 
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standard ET system-bath model Hamiltonian reads 

Ht = ha\a){a\ + {h + E°)\b){b\+Vi\a){b\ + \b){a\). (1) 

Here, E° denotes the reaction endothermicity, V the 
transfer coupling matrix element, and ha (hb) the sol- 
vent Hamiltonian for the ET system in the donor (ac- 
ceptor) state. The system is initially in the donor \a) 
site, with the solvent (bath) equilibrium density matrix 
Pa'^ oc e"''"/^'^^-^) at the specified temperature T. 

The ET system of eqlT] can be treated as a spin-boson 
problem in the context of quantum dissipation. We 
have recently constructed a general theory of quantum 
dissipation)2^i22ii^ which results in an analytical solution 
to the ET dynamics in a Debye (solvent) dissipatiou i^^i^^ 
It is noticed that the exact construction leads always to 
a generalized rate equation>^iS£i^ 



Pa{t) - - / drkit - T)PaiT) + / drk' (t ~ T)PbiT) . (2) 

Jo Jo 

Here, k{t) and k'{t) denote the forward and backward 
rate memory kernels. On the other hand, one often finds 
in practice that the simple kinetic theory, with a Marko- 
vian rate constant description, can well describe the ob- 
served rate process. The resulting reactant population 
Pa(t) decays exponentially toward its equilibrium value 
at a given temperature. Note that formally the rate con- 
stant is just the integrated rate kernel over time. Is there 



any quantitative justification for the Markovian rate pro- 
cesses being often observed experimentally? In this work, 
we try to address this issue, along with the quantum sol- 
vation effect and the adiabatic versus nonadiabatic na- 
ture of ET rate processes. 



B. Background 

Previous work on ET theory focuses mainly on the rate 
constant description. Consider Fig.[Tl the schematics of 
the simple donor-acceptor ET system of eq[TJ Here, the 
relevant (macroscopic) solvent potential surfaces Va and 
Vb are plotted as the functions of solvation coordinate 
U = hb — ha- The solvation energy for the ET from the 
donor \a) to the acceptor \b) site is given hy^^i^ 



A^tre {Up'a^)^{U). 



(3) 



The second identity here defines the notation (• • •) for 
the bath ensemble average, where the trace runs over all 
the solvent (bath) degrees of freedom. At the crossing 
{U + E° = 0) point, Va = Vb = {E° + A)V(4A). It is 
the celebrated Marcus' ET reaction barrier heightfii^i^ 
Thus, Fig.[T] also summarizes the Marcus' nonadiabatic 
rate expression. 



v^/h 



^JXk-BT/TT 



exp 



{E° + Xf 



dAfceT 



(4) 



This rate (constant) can be derived readily from eq[Tl 
by using the classical or static Franck-Condon approxi- 
mation, followed by the classical fluctuation-dissipation 
theorem, {U"^) — ([/)^ — 2k^T\. It does not consider the 
dynamic solvation effect, which is characterized by the 
relaxation time and also associated with the viscosity 
of the solvent. li>i2>i^>Hii5a6 Such effect was first stud- 
ied by Kramers in his classical Fokker-Planck-equation 
approach to the rate theory of isomerization reaction.— 
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The resulting rate shows the celebrated turnover behav- 
ior that the rate has a maximum in an intermediate vis- 
cosity region i^^'^^ This clearly demonstrates the dual role 
of solvent on reaction rate. Note that in the Kramers' 
rate theory, the reaction system is treated as a Brownian 
particle moving on a single adiabatic double-well poten- 
tial surface. 

The dynamic solvation on nonadiabatic ET has been 
incorporated in the quantum extension of Marcus' the- 
ory, formulated via the Fermi golden rule.— li^i^Sii^ The 
solvation coordinator is now a time-dependent stochastic 
operator U{t) = Q^hat/h-jj^-ih^t/h ^ assumed to follow 
the Gaussian statistics. Thus, the dynamic solvation is 
completely characterized by the correlation function. 



C{t~T) = {[U{t)-\][U{r)-\]). (5) 



The Fermi-golden-rule formalism is valid for an arbitrary 
form of the solvation correlation function. But like the 
Marcus' theory, it is restricted in the second-order trans- 
fer coupling regime. The resulting rate does not show the 
Kramers' fall-off behavior— i2S that involves the barrier 
recrossing and thus depends on higher-order transfer cou- 
pling. The improved approach has been proposed, on the 
basis of the fourth-order perturbation theory, followed 
by certain resummation schemes.— 

The resulting rate does support a smooth transition be- 
tween nonadiabatic and adiabatic ET reaction, recover- 
ing properly the Kramers' turnover behavior . ^^'^^ 

In this work, we adopt the reduced density matrix 
dynamics approach that is closely related to our recent 
development of quantum dissipation theoryi^i^i^i With 
the aid of the analytical expression for the nonperturba- 
tive and non-Markovian ET ratCfS^"^ we will elaborate 
in detail the quantum solvation effects in both the adi- 
abatic and nonadiabatic reaction regimes, and further 
investigate their relations to the Markovian versus non- 
Markovian nature of ET reaction kinetics. 

The remainder of this paper is organized as follows. 
Section reviews the generalized kinetic rate theory, 
constructed readily via the reduced density matrix dy- 
namics. In Sec. mil we decompose the total rate into its 
adiabatic and nonadiabatic components, and also an- 
alyze the Markovian versus non-Markovian nature of 
population transfer dynamics. We propose to use the 
Kubo's motional narrowing function, originally used in 
the context of optical spectroscopy i^^i'^^ to character- 
ize the Markovian character of ET rate process. Sec- 
tion IIVI presents a full account of the effect of quan- 
tum versus classical solvation on some representative 
ET reaction systems. Section |V] analyses the popula- 
tion transfer dynamics, with the focus on their adia- 
batic/nonadiabatic and Markovian/non-Markovian char- 
acters. Finally, Sec. lVIl concludes the paper. 



II. NON-MARKOVIAN RATES: THEORY 
A. Generalized rate theory and dissipation 

Let us define the dynamic rates via the evolution of 
reduced system density matrix p. In the absence of time- 
dependent external field, the generalized quantum master 
equation reads, 

p{t) - -jLp(i) - / dTii{t - r)p(r). (6) 
Jo 

Cp = hr^[H,p\, with H = (Ht) being the reduced sys- 
tem Hamiltonian; Il{t — r) denotes the dissipation ker- 
nel. Note that eqlHlis formally exact; in fact, the involving 
nonpcrturbative Tl{t) can be generally expressed in terms 
of continued fraction formalismj ^'^i'^^ '^ 

To obtain the dynamic rates, we shall eliminate the 
coherent components pjk', j ^ k, from eq[Sl and retain 
the populations Pj[t) = Pjj{t) only. To do that, let us 
recast eq[5] in its Laplace frequency-domain resolution, 

sp{s)~p{Q) = ~-[iC + n{s)]p{s), (7) 

and then arrange it into the block matrix form for the 
population and coherence vectors, P = {pjj} and C = 
{pjk',j 7^ k}, with C{t = 0) = 0, respectively. We have 

sP{s) - P(0) = -Tpp(s)P(s) - Tpc(s)C(s), (8a) 

sC{s) = -Tcp(s)P(s) - Tac{s)C{s), (8b) 

with the transfer matrices, Tpp, Tpc, Tcp, and Tec, aris- 
ing from the corresponding rearrangement of iL + n(s). 
Eliminating the coherent component C leads to 

sP(s) - P(0) = K(s)P(s), (9a) 

or, equivalently, 

Pi(t) = f '^^'^Mt - r)Pk{T). (9b) 
k -^0 

The involving rate resolution matrix is obtained as 

/•oo 

K(s)= / dte-^*K(t) =Tpc(s+Tcc)"'Tcp-Tpp. (10) 

Its element Kjk{s) resolves the statc-to-state dynamic 
rate memory kernel Kjk (t) . The rate matrix satisfies the 
relation Kjk = 0, as inferred from the population con- 
servation, Tip = J2j Pj = constant. 

B. Generalized rate theory in the two-state system 

For the two-state ET system (cf. eq[T] and Fig.[T]) of 
present study, only a single dynamic rate equation is in- 
dependent. That is eq[21 which reads in the Laplace do- 
main as 

sPais) - Pa{t = 0) = -k{s)Pa{s) + k'{s)Pi{s). (11) 
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Here, fc(s) = L{k{t)} = -Kaa(s) and k'{s) = L{k'{t)} = 
Ki,a{s) are the forward and backward rate resolutions, re- 
spectively. The involving transfer matrices for the simple 
ET system of eq[T]have all been identified; cf. the eqs 32 
of Ref. i25|: 



with 



Tp 



0, Tpc = — 



-1 1 
1 -1 



Here, 



i{E° 



-iV/h z* + iV/h 
iv/h z - iV/h 



X)/h y* 

x + i{E° + X)/h 



H 



ba,ba-} V 



H 



ha,ah-! 



z = H, 



ba.bbt 



(12a) 



(12b) 



(12c) 



(13) 



and their complex conjugates H*^-, = lij'j^k'k arc the 
only nonzero elements of the H-tensor. Denote also 

ais) = s + xis) + {i/h){E° + X). (14) 

The final expressions for the rate resolutions are then^^ 



k{s) 



2\Vl 

ft2 



Re 



ajs) + y{s) 
HsW-\yisW 



(15a) 



and 



[a{s)+y{s)][l-thz*{s)/V] \ 
H.sW-\y{sW J 



The involving parameters, defined in c q|131 can in prin- 
ciple be evaluated in terms of continued fraction formal- 
ism of the nonperturbative dynamics of reduced density 
matrixi^^i^i^ In particular, the continued fraction ex- 
pressions of these parameters have been solved analyt- 
ically for the Debye solvent model, where the solvation 
correlation assumes a single exponential form.— 



III. NATURE OF RATE PROCESSES: 
ADIABATIC VERSUS NONADIABATIC AND 
MARKOVIAN VERSUS NON-MARKOVIAN 

A. Adiabatic— nonadiabatic rate decomposition 



9(t) 



Ue° + X)t+^ fdT CdT'Cir'). (16b) 
n Jo Jo 



The corresponding A:na(s = 0) is the 
well established quantum nonadiabatic rate 
expression^iiaaiii2a3ii4a5a6^a8a9^^^ The Mar- 
cus' expression, eqlH can be obtained via setting 
C{t) w C(0) = 2XkBT that amounts to the static (slow- 
modulation) approximation, followed by the classical 
fluctuation-dissipation theorem. 

Consider now the long-time regime, involving only the 
integrated rate kernel, i.e., the rate constant k = k{s — 
0) or k' = k'{s — 0). Hereafter, rate constant will be 
referred in short as rate if it causes no confusion. For the 
later use, denote also the reaction equilibrium constant. 



Pb{oo) k{s = 0) 



Pa (OO) k'is = 0) 



k 



(17) 



It is noticed that k < /c^a is expected to hold in gen- 
eral; see Fig. inland its comments followed, especially on 
the case of exception. This inequality may be inferred 
from the fact that the total k contains the backscatter- 
ing contributions, while the nonadiabatic /c^a does not. 
The decomposition of the total rate into its nonadiabatic 
and adiabatic components may therefore be introduced 



a* 



11.12.13.14.15.16 



k\ 



(18) 



As k and fc^A can be evaluated via eq llSal and eqll6( re- 
spectively, with s — 0, the above equation can in fact 
be considered as the working deflnition of adiabatic kj^. 
We shall show later that fc^ is relatively much insensitive 
(compared with /cna) to the transfer coupling strength 
V variable (cf. Fig. [5]). Therefore, it can be practically 
used to describe the adiabatic rate process that involves 
only the ground solvation surface via the diagonalization 
of eq[Tl The ratio /cna/^a can be considered as the adia- 
baticity parameter. The ET reaction assumes adiabatic 
when k^j^/kA 3> 1, and nonadiabatic when k^j^/kA <ti 1. 



Characterization of Markovian versus 
non-Markovian rate processes 



It is noticed that the ET in the short-time (t < h/V) 
regime is always nonadiabatic, or fc oc F^. This fact can 
be inferred from the observation that the backscattering 
events, responsible by the higher-order transfer coupling, 
are yet not to occur. In Appendix|Xl we solve eqs[8] and 
1121 in the weak transfer coupling limit, and obtain the 
nonadiabatic counterpart of ea llSI 



k„A{s) = ~^R*^ J exp[-st - g{t)], 



(16a) 



We now turn to the Markovian versus non-Markovian 
nature of the reaction. The general theory of rate 
presented in the previous section assumes always non- 
Markovian. On the other hand, the experimental obser- 
vations appear often Markovian. It is desirable to have a 
working criterion on the nature of ET kinetics. In contact 
with experiments, let us consider the scaled population. 



A(t) 



P,(0)-P,(oo) 



(19) 
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It does not depend on the state-index due to the identity 
of Pa{t) + Pi,{t) — 1. The simple (Markovian) kinetic rate 
equation, Po(0 = —kPa{t) + k'Pb{t), can be represented 
in terms of the scaled population as 



(20) 



The involving decay constant, w = k + k', could be 
measured readily, if the rate does behave like Marko- 
vian. The forward and backward reaction rate con- 
stants can then be evaluated as fc = wKcq/{l + K^q) 
and k' — w/{l + K^q), respectively. 

A non-Markovian rate process can be described by the 
kinetic rate memory kernel w{t — r) via 



A{t) 



dTw{t - t)A{t). 



(21) 



In other words, one can deduce the rate kernel from the 
population evolution as w{t) = L~^{[1 — sA{s)]/A{s)}. 
Here A(s) denotes the Laplace transform of A{t), while 
L~^{-} the inverse Laplace transform, i.e. 



w{s) = L{w{t)} = dte 
Jo 



^-^(i)^i_f^. (22) 
A(.) 



In AppendixiBl we present the explicit expression of w{s) 
in terms of k{s) and fc'(s); see ea lB7[ together with some 
useful identities in relation to the Laplace transform. Un- 
like w(s), the population evolution alone is in general not 
sufficient to determine both k(s) and k'{s). 

Equivalent to the kinetic rate kernel w{t) defined in 
ea[2Tl one may also define the time-local rate W{t) 
via A{t) = -W{t)A{t). Note that while the kernel 
w(t) is always well behaved, its time-local counterpart, 
W{t) = — A(t)/A(i), may diverge at certain time, say 
t', due to the possibility of A(t') = in an under- 
damped non-Markovian rate process. At time t' , the 
effective kinetics reduces to the zero-order rate process, 
since W{t')A{t') remains finite. Moreover, eq[2l1 implies 
also that W{t = 0) = and W{t oo) ^ w{s = 0) = wq- 
The latter is nothing but the fact that long-time regime 
is Markovian. 

To determine the short-time behavior for the popula- 
tion dynamics, we notice the identity. 



k{t^O) ^ k'{t = Q) = 2V^/h^ 



(23) 



This can be obtained directly by considering the fact that 
the short-time behavior is identical to ki^A{t = 0) (cf. 
eg fTB)) . regardless whether the reaction is nonadiabatic 
or not. Together with eq |B8[ and setting Pa(0) = 1, we 
obtain 



w = w(t^O) = 



(24) 



oq 



It determines the short-time behavior of the scaled pop- 
ulation dynamics, as A(0) — and A(0) = —w that are 
implied in ea l21l 



The short- and long-time behavior described above can 
be summarized as 



A(t->0) ^e-'^*'/^, A(t->oo) 



(25) 



Besides the above asymptotic behaviors, A{t) is in gen- 
eral also influenced by the coherent motion or quantum 
beat, as long as the rate process goes beyond the second- 
order in the transfer coupling V . As the asymptotic be- 
haviors are concerned, it may suggest to use the Kubo's 
line shape function,— 



A^it) = exp 



-2/ -wt/w 



1 -f wt/w] 



(26) 



to analyze the nature of population transfer dynamics. 
The involving Kubo's Markovianicity parameter is (cf. 
eqs[2land[25l) 



2K, 



cq 



1+K, 



cq 



1/2 



V 



(27) 



The rate process assumes Markovian or non-Markovian, 
when K > 1 or K < 1, respectively. 



IV. EFFECTS OF QUANTUM SOLVATION: 
NUMERICAL RESULTS AND DISCUSSIONS 

A. Debye solvent model and general remarks 

For the numerical demonstrations presented hereafter, 
we consider the ET system in eqU] or Fig.[l] with the 
solvation correlation function in eq[5] being characterized 
by 



C(t) =77CXp(-t/TL). 



(28) 



This is the Debye solvent model, with the longitudinal re- 
laxation time tl and the pre-exponential parameter— 



r]^Xi2kBT-ih/n^). 



(29) 



For the above ET system, the analytical expressions for 
the reduced density matrix p{t) and the involving dy- 
namics rate functions, k{s) and fc'(s), have all been con- 
structed in Ref. 25. Note that eq[28l satisfies a semiclassi- 
cal fluctuation-dissipation theorem, valid when the tem- 
perature is comparable with or higher than the system's 
transition energy .^^'^^ This is the only approximation in- 
volved in this paper. 

To elucidate the quantum solvation effect on ET, es- 
pecially as the reaction mechanism is concerned, the rate 
in the classical solvation limit will also be evaluated as a 
reference. The classical solvation correlation function as- 
sumes real; i.e., i] ijd = 2AA;bT when ksT/h » t^^ . It 
follows \r] - 77ci|/??ci = 0.5rther/rL, with rthor = ^/{k^T) 
denoting the thermal time (rther = 26 fs for T — 298 
K). The quantum solvation effect can be significant when 

Tl < Tthor- 
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It is also noticed that the solvation longitudinal relax- 
ation time tl is proportional to the solvent viscosity Ji^ii^ 
In this sense, the k versus tl behavior can be referred 
as the rate-viscosity character. This connection suggests 
the well-established Kramers' picture of the solvation ef- 
fect on chemical reactio n^^'^^ be exploited in the follow- 
ing to elucidate the underlying ET reaction mechanism. 



B. Quantum vs. classical solvation effects 

As the mechanism is concerned, the effect of quantum 
versus classical solvation is expected to be most distinct 
in the following two scenarios (cf. Fig.[T]). One is the sym- 
metric ET system {E° ~ 0) in which the Fermi resonance 
enhanced quantum tunneling is anticipated. Another is 
the classical barrierless system {E° -I- A = 0) where the 
Marcus' inversion takes place. Figure [2] depicts the eval- 
uated rate-viscosity characteristics for these two scenar- 
ios of ET system, at two values of transfer coupling: 
V = IkJ/mol (upper-panels) and V = O.OlkJ/mol 
(lower-panels). The solvation energy A = 3 kJ/mol and 
temperature T = 298 K. 

In the high-viscosity (tl > Tthor) regime, the differ- 
ence between the quantum (solid-curve) and the classical 
(dashed-curve) in each panel diminishes (at the qualita- 
tive level). This observation is consistent with the phys- 
ical picture that the high viscosity (or slow motion) im- 
plies a large effective mass and thus leads to the classical 
solvation limit. Observed here is also the Kramers' fall-off 
behavior^^-^^ for the adiabatic rates in the upper panels. 
This is the diffusion limit; the higher the solvent viscosity 
is, the more backscattering (or barrier re-crossing in the 
classical sense) events will be. For the nonadiabatic pro- 
cesses in the lower panels, the backscattering effects are 
quenched. This accounts for the plateau, observed in the 
lower-panel (c) or (d), at which the Marcus' nonadiabatic 
ET regime is reached, and the rate becomes independent 
of viscosity. 

In the low-viscosity (tl < Tthcr) regime, the difference 
between quantum and classical solvation is at the mech- 
anism level for the two specified types of ET systems. 
For the symmetric {E° = 0) system [the left-panel (a) 
or (c) of Fig. [2], the quantum rates (solid-curves) are ap- 
parently Fermi resonance-assisted tunneling dominated 
processes, while the classical rates (dash-curves) are bar- 
rier crossing events. In contrast to that k 3> fcd in the 
Th < TthcT regime for the symmetric {E° = 0) system, 
observed is the opposite result of fc <C fed for the Marcus' 
barrierless {E° -I- A = 0) system [the right-panel (b) or 
(d) of Fig. [2]. In the latter case, the classical rate does 
behave as a barrierless reaction, but the quantum rate 
exhibits the Kramers' barrier-crossing characteristics as 
a function of viscosity. This may indicate that for the 
ET in the classical barrierless system there is an effective 
viscosity-dependent barrier that vanishes as tl increases. 
In other words, the barrier of (£'° + A)^/(4A), as depicted 
in the ET schematics Fig.[Tl is the static picture that 



assumes the solvation potential being fixed. This pic- 
ture is valid in the high-viscosity (or slow-modulation) 
regime, but no longer true in the low-viscosity (or fast- 
modulation) regime. In general, the effective solvation 
potential for ET reaction is solvent viscosity dependent. 

To confirm the above observed ET mechanism-related 
features, the values of A and T are varied, and the results 
for V = IkJ/mol are summarized in Fig. [31 Here, k/kd 
as functions of TL/rthcr are depicted. This figure verifies 
that TL/rthcr does serve a proper measure for the nature 
of solvation. The reaction mechanism turnover occurs at 
Th — Tthor, for the symmetric {E° = 0) case, and classical 
barrierless {E° -I- A = 0) ET systems. In the former case, 
it is changed from the tunneling to barrier-crossing, while 
in the latter case from the barrier-crossing to barrierless 
ET rate process. 

Let us now consider the Marcus-type inversion behav- 
iors. Figure m presents the Marcus-plots, the logarithmic 
rate log k versus reaction endothermicity E° , in relation 
to Fig. [2] Three values of TL/rthcr are chosen as 0.1 (solid- 
curves), 1 (dot-curves), and 10 (dash-curves), to repre- 
sent the low-, intermediate-, and high-viscosity regimes, 
respectively. 

Note that E° — —X does represent the classical barri- 
erless ET systems in all cases in study. All the classical 
rates [right-panels, Fig.|4l^b) and (d)] have inversions oc- 
curring at E° — ~X. Moreover, all of them are symmetric 
about E° = —A. The Marcus' parabolic character is re- 
covered in the high-viscosity, nonadiabatic and classical 
limit, and is practically the dashed curve in Fig.^d). 

The quantum rates depicted in the left-panels [Fig.[3]^a) 
and (c)] also show inversion behavior. However, the ob- 
served inversion region depends sensitively on the solvent 
viscosity, and shows also an asymmetric character that 
will be elaborated soon. Apparently, the quantum rate 
inversion region is closely related to the interplay between 
the barrier-crossing and tunneling processes. 

Let us start with Fig.[4Ua). Consider first the high- 
viscosity regime, in which, according to the analysis made 
earlier for Fig. [21 the ET process is qualitatively the same 
as its classical counterpart. Consequently the dash-curve 
in Fig. Ufa) has the inversion region around the classical 
barrierless position at E° = —A. Consider now the low- 
viscosity regime, in which, again, according to the anal- 
ysis earlier there is always a nonzero barrier for the ET 
reaction, covering over the entire range of E° including 
the value of E° — —A. This explains the inversion be- 
havior of the solid curve of Fig.[3I^a) that is peaked at the 
resonant position of E° — 0. As the viscosity increases, 
the inversion region smoothly shifts from the resonant 
peak position E° = to the classical barrierless position 
ofE° = -A. 

To explain the asymmetric property of the quantum 
inversion behavior as depicted in the left panels of Fig.[31 
recall that k{—E°) « k'{E°), the backward reaction rate, 
and k{E°) < k'{E°) for an endothermic {E° > 0) reac- 
tion. This leads immediately to the asymmetric property 
of the solid-curve in Fig.^a) or (c), in which the blue 
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(endothermic) wing falls off faster than its red (exother- 
mic) wing. This asymmetry decreases as the viscosity 
increases, since the high viscosity regime behaves classi- 
cally. 



V. THE ADIABATIC AND MARKOVIAN 
CHARACTERS OF THE REACTION 

We are now in the position to demonstrate (for the 
cases of quantum solvation only) the issues in relation to 
the nature of ET rate process, as addressed in Sec. lIIII 
Let us start with the adiabatic/nonadiabatic character 
depicted in Fig.O Here, each individual k (thin)-curves, 
as function of tl, is decomposed into its adiabatic fc^ and 
nonadiabatic /cna components (e allSp . These two compo- 
nents are given in the upper-panel [Fig.[5Ka) or (b)] and 
the lower-panel [Fig.EKc) or (d)], respectively. The left- 
panels are for the E° = system [Fig.[SJa) and (c)], while 
the right-panels are for the + A = system [Fig.tSfb) 
and (d)]. In each panel, four values of transfer coupling 
are used: y/(kJ/mol) = 0.25 (dotted), 0.5 (dashed), 1 
(solid), and 2 (dash-dotted). 

It is observed that in the low- viscosity (tl < Tthor) 
regime, the reaction is nonadiabatic; see Fig.[5I^c) and 
(d). This observation may be understood as follows. 
In the low-viscosity regime, the solvent fluctuates fast 
and stabilizes the ET system in the acceptor state before 
backscattering taking place. As results, the reaction is 
nonadiabatic in the low-viscosity regime, at least for the 
range of transfer coupling strength considered here. The 
above picture is also consistent with the observation that 
in the diffusion (tl ^ Tthor) limit, where k oc 1/tl, the 
reaction assumes an adiabatic rate process; see Fig.[5lja) 
and (b). 

Figure [S] presents the adiabaticity parameters, /cna/^a, 
as function of reaction endothermicity E° , at V = 1 
kJ/mol, with the three specified values of TL/rthcr = 0.1, 
1, and 10. The inverse of these values are also used indi- 
vidually to scale the corresponding adiabaticity curves, 
as depicted in Fig.[6l The adiabaticity curve in high- 
viscosity (rL/rthor = 10) regime looks all normal, being 
of the minimum about where the rate maximum is; cf. 
the dash-curve in Fig.^a). The faster the ET passage, 
the less adiabatic (or more surface-hopping) the reaction 
would be. This feature remains largely unchanged in 
the intermediate- viscosity (rL/rther = 1) case, except for 
that the adiabaticity minimum now is slightly negative, 
or fc > fefjA- It contradicts with the argument made ear- 
lier for e q|18l The abnormality here may be accounted 
for, at least partially, by the associated tunneling rate 
process, as elaborated below. 

The abnormality in the adiabaticity parameter is most 
striking around the symmetric {E° « 0) system in the 
low-viscosity (rL/rthor = 0.1) regime. This is tunnel- 
ing dominant scenario, as discussed in Sec. lIVI Con- 
sidered here is also the strong transfer coupling case of 
V = IkJ/mol. Thus, the observed abnormality is most 



likely caused by the coherent ET reaction. The argument 
for k < /cnaj due to the backscattering- induced total rate 
reduction, may no longer be valid. Moreover, we will 
see soon that the ET reaction in the observed abnormal 
region is non-Mar kovian, leading to the rate constant de- 
scription inadequate at all; see Fig. [TO] and the comments 
there. 

We now turn to the Markovian/non-Markovian nature 
of ET reaction. The key quantity here is the Marko- 
vianicity parameter k fea l27p . It involves the reaction 
rate k and the equilibrium constant K^q. The detailed 
knowledge on how i^oq depends on the Debye solvent pa- 
rameters can be found in Ref. |2^ For the rate in the 
low-viscosity regime, k « fc^A, as depicted in the lower- 
panels of Fig.O resulting in k x 2V'^V/{E°'^ + F^), with 
F = 2\k-QTTi^ for the Debye solvent model in study; cf. 
eqlHwith eq[28l Thus, 

K « Xk^TTi^/V; when E° and tl < nher- (30) 

This equation can be used to estimate the Markovianicity 
parameter for a symmetric ET system {E° = 0, implying 
also = 1) in the low-viscosity regime. 

Figure [7] shows the evaluated Markovianicity k fea [77|) 
as function of rL/rthor, for the same ET systems as Fig. [2] 
(without the classical solvation parts). For E° + X = 
(lower-panel), the ET rate process behaves Mar kovian 
(k > 1), even for the strong transfer coupling {V = 
IkJ/mol) case. Apparently, the non-Markovian (k < 1) 
rate process is most likely to occur in the Fermi-resonance 
tunneling regime, where E° = and rL < rthor, with the 
approximated expression of k given in ea l30l The above 
comments are further confirmed by Fig.[Hl in which the 
Markovianicity k is plotted as function of E° , at the three 
representing values of rL/rthor- 

Figure [5] depicts the scaled population evolution, A{t) 
(eqlHD, for the ET systems of = (three left-panels) 
and E° + X = (three right-panels). The upper [(a) and 
(b)], middle [(c) and (d)], and lower [(e) and (f)] pan- 
els are of the specified values of viscosity, rL/rthor = 0.1, 
1 and 10, respectively. The transfer coupling strength 
is V = I kJ/mol. The weak transfer couphng {V = 
O.OlkJ/mol) counterparts, as depicted in the inserts of 
individual panels, are shown all Markovian, with k 3> 1 
for their values of Markovianicity; cf. Fig. [7] and Fig.[51 
Included for comparison in each panel are also the Kubo's 
AK(t) fea l^S)) and Markovian AMar = exp(— wi) counter- 
parts, where w = k + k' . As the envelop of population 
evolution is concerned, the significant non-Markovian na- 
ture is only observed in Fig.[HKa), with the relevant part 
is enlarged in Fig.[TOl This is a symmetric ET system in 
the strong transfer {V = IkJ/mol) coupling and low- 
viscosity (rL/rthor = 0.1) regime. The corresponding 
Markovianicity value (k = 0.3) is found to agree well with 
the aforementioned approximate expression of e ql30l 
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VI. CONCLUDING REMARKS 

The quantum solvation, adiabatic/nonadiabatic and 
Markovian/non-Markovian characters are important is- 
sues in understanding chemical reaction, including ET 
in solution. Here, we have revisited these issues in a uni- 
fied and transparent manner, with the aid of Debye sol- 
vent model (eq[28l) that supports an analytical solution 
without additional approximations!^ The physical pic- 
ture discussed in this work is however rather general. 

We have presented a full account for the effect of quan- 
tum solvation on the ET rate process fSec. lIV|) . Not 
just can it change a barrier-crossing event to tunnelling, 
the quantum nature of solvent can also lead a classical 
barrierless reaction to an effective barrier-crossing rate 
process. The resulting rate may differ from its classical 
counterpart by order of magnitude. The quantum solva- 
tion is found to be distinctly important in low-viscosity 
(fast-modulation) solvents. For a realistic solvent that 
consists of multiple correlation time scales, only the slow- 
modulation solvent modes can be treated classically. 

The adiabatic-nonadiabatic decomposition of rate 
feg fTS)) is practically useful, provided that the total rate 
k can be experimentally measured and the nonadia- 
batic rate fc^A can be readily evaluated via ea ll6l The 
adiabaticity parameter /cna/^a — ^na/^ ^1, as in- 
ferred from e gflSl can then be used to discuss the adi- 
abatic/nonadiabatic nature of reaction. Interestingly, a 
negative adiabaticity may indicate there is a certain de- 
gree of quantum tunneling taking place; cf. Fig. [5] and its 
comments. 

We have also proposed to use the Markovianicity 
parameter k (ea \T7]i . based on the Kubo's motional 
narrowing function, for analyzing the Markovian/non- 
Markovian nature of electron transfer rate process. Note 
the solvent relaxation time scale is typically in the order 
of picosecond, while Tthcr for room temperature is 26 fs. 
This amounts to the high-viscosity or slow-modulation 
regime of present studies on ET. The resulting Marko- 
vianicity parameter, as depicted in Fig. [7] and Fig. [SI is 
typically of k > 1 for the aforementioned typical cases. 
This may account for why most experimental observa- 
tions do support the Markovian [k > 1) rate constant 
description. 

We have pointed out that a non-Markovian rate pro- 
cess is most likely to occur in the symmetric ET system 
at the fast-modulation regime. It is just the opposite to 
the spectroscopic case. According to the motional nar- 
rowing picture, the fast modulation leads to a Markovian 
spectroscopic process. "^^'^^ The above seemingly counter- 
intuitive phenomenon in relation to the nature of rate 
process may be understood as follows. First of all, the 
motional narrowing picture is applicable to the spectrum 
of rate kernel, rather than the population evolution it- 
self. The narrower the rate kernel spectrum is, the less 
Markovian of rate process would be. However, this is not 
the complete picture. The peak position of the rate ker- 
nel spectrum, in relation to where rate constant is eval- 



uated, should also be considered. It shifts from the clas- 
sical Marcus' inversion position at E° = —X in the slow- 
modulation limit, to the quantum resonant tunneling at 
E° = in the fast-modulation (or motional-narrowing) 
regime, see Fig.|4] (left-panels). Equation ([30)) that is 
achieved at E° — can be considered as the lower bound 
of the Markovianicity k for the ET rate process, in the 
fast modulation regime. In this regime, the population 
transfer may also exhibit the quantum beat feature that 
is non-Markovian in a strict sense. We shall investigate 
these complex cases elsewhere. 
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APPENDIX A: PERTURBATIVE RATE 
KERNELS 

In this appendix, we shall treat the nonadiabatic rate 
problem, on the basis of the standard perturbation the- 
ory on the reduced density matrix p(t), assuming the 
transfer coupling last term of eq[TJ is weak. At the ini- 
tial time t = 0, the total composite density operator is 
Pt(0) — p'^\a){a\. Consider eq[8|for the coherence com- 
ponents, which for the two- level ET system reads explic- 
itly as (setting h = I in this appendix) 

apba = -iViPa - Pb) - Zpb - Vpla- (Al) 

Here a — s + i{E° + X) + x, and 

X = Uba,ba, y = Hba.afc, Z = Ilba,bb- (A2) 

Combining the initial condition Pa(t ^ 0) = 1 and the 
perturbative action of V together result immediately in 
Pi°^ = 1/s and 

^ p{2k+l} ^ p{2k} ^ ^{2fc+l} ^ ^^{0} ^ y{2k+l) ^ ^{2k}^ 

Therefore, the lowest order in ea lAll reads 

aW^i} = -^V/s. (A3) 

On the other hand, the standard first-order perturbative 
expression is 

p\^^{t)^-iV fdreM-gir)]- (A4) 
Jo 

Here 

exphg(i)] = ^exp+ [-i j^dT [E° + C/(t)]|^ . (A5) 
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Using the second-cumulant expansion expression, which 
is exact for the Gaussian solvation process, results in 

g{t) = i{E° + X)t+ f dT f dr'C(T'). (A6) 

We have then 

spllHs) = ^iVL{eM-9im ^ -iVJis). (A7) 

Together with eq lA3| we obtain 

= 1/J(s). (A8) 

Together with y^'^^ — 0, ea llSal to the lowest order 
reads k^A{s) = {2V^ / h?)ReJ {s) , where J(s) denotes the 
Laplace transform of exp[— g(i)] (cf. eq [A7|) . We obtain 
therefore e afTCl 



APPENDIX B: SOME USEFUL RELATIONS FOR 
RATES 

Let us first present some basic relations in connection 
to the Laplace transform, defined for s > as 



(B3) 



f{s)^L{f{t)}^ / dte-^'f{t). 



It satisfies the boundary condition, L{-}\ 



dTh{t~r)h{T)\^h{s)h{s), 



(Bl) 
0, and 
(B2) 



Using eq lB3| together with the identities of ~ 

and L{/(i)}|s=o ~ /(oo) ^ /(O)? '^6 obtain immediately 



/(O) = lim [sf{s)l /(^) = lim[s/(s)]. (B4) 

s — >oo s — >0 



We are now in the position to derive some useful rela- 
tions between the non-Markovian rate variables appear- 
ing m SecHIl From eqEland Pa{t) + Pb(i) = 1, we have 



Pa{3) 



Pa{0) + k'{s)/s 



s + k{s) + fc'(s) ■ 
From eqs[TH] and [221 we have 

Pais)- Pa{00)/S 1 



PaiO)-Pa{oo) S + W{s) 

The above two equations lead to 

s[Pa{0)k{s)-Pb{Q)k'{s)] 



(B5) 



(B6) 



w{s) 



s[Pa{0) - PaM] - PaMfc(s) + Pb(oo)fc'(s) ' 

(B7) 

Together with the first identity of eq lB4[ we have 



PaiO)kit^O)-Pbm'it = 0) 

= = P.(0)-P.(oo) ■ 

The above relation will be used in deriving eq |24l 
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FIG. 1: Schematics of solvent potentials Va and Vb for the 
ET system in the donor and acceptor states, respectively, as 
fimctions of the solvation coordinate U = hb — ha = Vb — Va — 
E°, with E° being the ET endothermicity and A = {U) the 
solvation energy. The classical barrierless system is that of 
S° + A = 0. 



FIG. 2: Electron transfer rate k (solid-curves), as function 
of solvent longitudinal relaxation time tl , with A = 3 k J/ mol 
at T = 298 K. Left-panels (a) and (c): E° = 0; right-panels 
(b) and (d): E° + \ = 0; upper-panels (a) and (b): V = I 
kj/mol; lower-panels (c) and (d): V = 0.01 kj/mol. Included 
in each panel is also the classical solvation counterpart fed 
(dash-curves). Note that Tther = fi/{kBT) — 10~^'®ps. 



FIG. 9: The scaled population A{t) fe all9|) evolution, eval- 
uated at TL/rthor =0.1 (upper-panels), 1 (middle-panels), 
and 10 (lower-panels); Left-panels: E° = 0; right-panels: 
E° + \ ^ 0. Here, V = IkJ/mol, A = 3kJ/mol, and 
T — 298 K. Included in each panel are also the corresponding 
Kubo's Ajiit) (e Ql26l dashed) and Markovian AMar(i) (dot- 
ted). The insert in each panel is the V = O.OlkJ/mol coun- 
terpart, where all these three curves are identical. 



FIG. 3: The ratio of quantum versus classical rates, k/kd, as 
function of TL/rthcr, with the transfer coupling V — 1 kJ/mol, 
at various specified values of A and T. 



FIG. 4: The rate k (left-panels) and the classical counter- 
part feci (right-panels), as function of E° , with A = 3 kJ/mol 
at T = 298 K; upper-panels (a) and (b): V = 1 kJ/mol; 
lower-panels (c) and (d): V — 0.01 kj/mol. The three val- 
ues of relative relaxation time scale, TL/rthor ~ 0.1 (solid), 
1 (dotted), and 10 (dashed) are chosen to represent the low, 
intermediate, and high- viscosity regimes, respectively. 



FIG. 5: The decomposition of rate k (thin-curves), following 
e allSI into the adiabatic fcA (upper-panels) and nonadiabatic 
^NA (lower-panels) components, and plotted as functions of 
TL/rthor, with A — 3kJ/mol and T = 298K. Left-panels (a) 
and (c): E° = 0; right-panels (b) and (d): £;° -|- A = 0. 
Each panel involves four values of transfer coupling strength: 
l//(kJ/mol) = 0.25 (dotted), 0.5 (dashed), 1 (solid), and 2 
(dash-dotted), respectively. 



FIG. 6: The adiabaticity /cna/^a as function of E° , for 
Th/ Tther = 10, 1 and 0.1. The inverse of TL/rther is also used 
to the scale individual adiabaticity curve. V = IkJ/mol, 
A = 3kJ/mol, and T = 298K. 



FIG. 10: The amplified portion of Fig.j^J^a). The scaled pop- 
ulation evolution A(t) (solid), Kubo's AK(i) (dashed), and 
Markovian AMar(^) (dotted). 



FIG. 7: The Markovianicity parameter k fe al27|l . as function 
of rL/rthor, at = 1 kJ/mol (solid) and 0.01 kJ/mol (dashed), 
with A = 3kJ/mol and T = 298 K, for (a) E° = and (b) 
£'° + A = 0. 



FIG. 8: The Markovianicity parameter k fe al27|) as the func- 
tion of E°, with A = 3 kJ/mol and T = 298 K, at rL/rthor = 
0.1 (solid), 1 (dotted) and 10 (dashed): (a) V = 1 kJ/mol 
and (b) V = 0.01 kj/mol. 
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